SELF-REGULATING GALAXY FORMATION 
I. HII DISK AND LYMAN ALPHA PRESSURE* 

Donald P* Cox^ 

Space Physics and Astronomy Department 
Rice University 


*Formerly: An Eddington-Like Limit to the Surface Density of a 

Galactic Disk 

^On leave from the Department of Physics, University of Wisconsin, 
Madison 



2 


Abstract 

Assuming a simple but physically based prototype for behavior 
of interstellar material during formation of a disk galaxy, 
coupled with the lowest order description of infall, a scenario is 
developed for self-regulated disk galaxy formation. Radiation 
pressure, particularly that of Lyman a (from flourescence 
conversion of Lyman continuum), is an essential component, 
maintaining an inflated disk and stopping infall when only a small 
fraction of the overall perturbation has joined the disk. 

The resulting "galaxies” consist of a two dimensional family 
whose typical scales and surface density are expressable in terms 
of fundamental constants. The model leads naturally to galaxies 
with a rich circuragalactic environment and flat rotation curves 
(but is weak in its analysis of the subsequent evolution of halo 
material). 

The galactic family has a natural lower limit on radius 

around 150 pc, at which point the systems are spheroidal and 

resemble giant HII regions. The largest radius is roughly 20 kpc, 

depending somewhat on the epoch of formation. The model seems to 

apply most naturally to large spirals with radii from 4 to 20 kpc 
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and average column (or surface) densities from 100 to 500 pc 

The formation scenario for a Milky-Way-like galaxy would be 
roughly as follows. The material eventually incorporated into the 

Q 

disk reached a maximum extension of 45 kpc at t ^ 5 x 10 years 


ABB (= After Big Bang). It was sufficiently sheared that the 

highest velocity material was only marginally bound. During 

recollapse, angular momentum was azirauthally shared, dictating a 

final radius smaller by a factor of four. Disk formation 

9 

commenced at t ^ 10 years, achieving an average surface density 
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of 300 pc only 5 x 10 years later. At this point, star 

formation is so rapid in the disk that Lyman a pressure dominates 
the disk structure as well as supporting the residual halo against 
further infall. The Lyman a supported halo moves toward 
hydrostatic equilibrium with M(r) * r. The fragmentation of the 
halo into bits (many possibly substellar) under the influence of 
Lyman a (or supernovae) would leave a galactic system with some 8% 
of its mass in the disk, the rest in a massive halo, possibly 
including the Magellanic clouds. 

The formation of such a system would be an extremely dramatic 
event. The peak luminosity would be of order 10^^ erg s~^, 

Lyman a accounting for 10 to 20%. It would last roughly 4 x 10^ 
years during which much of the disk material enters the first 
generation of stars and while several percent of the galactic 
hydrogen is burned to helium or beyond by stars with main sequence 
lifetimes less than that duration. The scale height of the ISM 
during formation would be about 150 pc but the lack of definiton 
of the mean plane would give the first generation of stars an 


inflated final distribution 



I. Introduction 


Any attempt to explore the formation epoch of spiral galaxies 
is naturally frustrated on four counts. We need to understand the 
universe well enough to predict the distribution of infalling mass 
in angular momentum and time. We need to understand the 
structure, phases, and mechanisms of the nascent interstellar 
medium (NISM). We need to be able to determine the star formation 
rate and zero-age initial mass function (ZAIMF) as functions of 
themselves and their interactions with the NISM. Finally we need 
to understand the feedback effects of the forming system on its 
own continued infall rate (or that of its neighboring systems). 

At the present time none of these four elements seem to be 
amenable to a confidence-inspiring analysis. 

Consider in particular the problem of determining the rate of 
star formation. Its solution cannot be divorced from 
understanding mass flow through the phases of the ISM. Mass 
filters from low density to high. It piles up at stable phases, 
accumulating until it achieves an adequate transfer rate to the 
next denser phase. Calculating the rate of mass transfer between 
phases and from the densest phase into stars involves 
understanding the heating, cooling, and environmental influences 
for each phase. In short, one must understand bottlenecks in the 
mass transfer. 

These bottlenecks are presently not well understood even in 



the solar neighborhood. We know that there are at least two 
stable phases, diffuse and molecular clouds, and something about 
each. It may be that diffuse clouds are heated by photoelectric 
ejection off dust grains by starlight. It may be that they are 
held together by the pressure of a hot low density matrix 
maintained by supernova explosions. But what governs the rate of 
cloud destabilization that diverts their mass to the denser 
molecular cloud phase? Do we even know for sure whether the 
passage is fairly direct or whether the diffuse clouds dissipate 
and regroup many times before entering the molecular cloud phase? 

Molecular clouds are largely opaque due to dust grains, and 
self gravitating. They may be heated in part by cosmic rays and 
by dissipation of turbulence of as yet unknown origin. They may 
be supported in part by magnetic fields, the latter kept tied to 
the material by ionization due also to cosmic rays. Electrons 
from elements with first ionization potentials below that of 
hydrogen may at times also be important. 

In contrast to this partial understanding, we almost totally 
lack knowledge about the ISM of a galaxy less than 10^ years old, 
the NISM. Such a system may offer comparable complexity with 
different forms, or it may offer a remarkable simplicity. It is 
quite possible that for a while there are no stable diffuse or 
molecular clouds. The heavy element abundance is very low, 
lowering dust content, cooling cfoef f icients , heating rates, and 



opacities. The galactic magnetic field may not yet have arisen. 
Cosmic ray acceleration or trapping may be inefficient. If so, 
what can we then infer about the NISM state and its star formation 
rate? 

One answer to the above question was proposed by Cox (1983, 
Paper I). It was suggested that the densest stable phase 
available to the NISM was ionized hydrogen at 10^K, with heating 
and ionization provided by stellar UV. This simple possibility 
leads to a trivial equation of state for the medium and a 
prescription on the formation rate of massive stars. The basic 
idea is that any gas which is not kept photoionized cools rapidly 
and forms stars, some of which are hot enough to contribute 
ionizing photons. The critical star formation rate is that which 
just manages to keep the residual medium fully ionized. 

In order to calculate the critical star formation rate, one 
needs to know only the efficiency of stars in producing ionizing 
photons. For timescales longer than the lifetime of the most 
massive stars, about 3 x 10^ years, an instantaneous recycling 
approximation is appropriate: UV photons are returned immediately 

from any new generation of stars, making the UV production rate 
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per unit volume where p* is the contemporary density per 

second going into stars and <j> is the number of ionizing photons 
returned per gram of star formation. 


* Summarizing from Paper I, we can write 
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where ra^ is the proton mass, 7 MeV is the energy liberated per 
proton burned, <hv> is the average ionizing photon energy, is 
the fraction of the hydrogen which gets burned after entering a 
massive star, is the fraction of the luminosity beyond the 
Lyman limit, and f^ is the fraction of the total mass in star 
formation which enters shortlived massive stars. Paper I 
estimated a consistent set of these parameters to be f^ = 0.5, f 2 
= 0.3, f^ = 0.06 I and <hv> « 20 eV, making m^<j) » 30001 photons 
per proton. (A slight correction due to the dilution by 
primordial helium will be introduced shortly.) The uncertain 
parameter I depends on the IMF. It is approximately 1 for the 
solar neighborhood IMF and increases with any weighting of the 
ZAIMF toward more massive stars. 

In what follows it will be assumed that the primordial ratio 

of helium to hydrogen nuclei was 0.1. The gas nuclei number 

density is n = n^ + - 1.1 whereas the gaseous mass density 

is p = 1.4 m = (1.4/1 ,l)ra n. It will also be assumed that 
g P M P 

both hydrogen and helium are singly ionized, making n 0 = n and 
that the recombination coefficient (to n > 2) a » 2 x 10 ^cm^s~* 
applies to both. (This coefficient is slightly lower than that 
for a present day HII region because the temperature is somewhat 



higher*) Thus the recombination rate per unit volume is an and 
the critical rate of proton flow into star formation is an /(m^), 
providing one ionization per recombination. Then we have 
(including helium) 
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For a volume average rate, the right hand side must be divided by 
the square of the filling factor f Q to get the local density from 
the average density, and then multiplied by the filling factor to 
count only the active volume. The net result is that the volume 
average star formation rate and gas density are related by 


= Ap, 


where 
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Two further points are made in Paper I. The first is that 
this rate probably represents an extreme upper limit to quasi- 
steady formation of stars. The inclusion of any other bottlenecks 
should provide a slower rate. The second is that when this rate 
is compared numerically to the rate which Larson (1969) found was 
needed to provide acceptable models for elliptical galaxies, the 
two are equal within uncertainties. 



The remainder of this paper takes the first steps at 
incorporating this specific NISM and star formation model into a 
scenario for the formation epoch of spiral galaxies. These steps, 
numbered by section, include: 

II. Evaluation of the structure, star formation timescale, 
and luminosity of a self-gravitating isothermal disk. 

III. Recognition of the importance of radiation pressure in 
disks of interest and a calculation of its magnitude, particularly 
that due to trapped Lyman a. Discussion of the nature of 

Lyman a pressure and its effects on the disk structure. 

IV. Discussion of the possibility that Lyman a pressure might 
seriously affect the infall rate, and determination of the 
properties of a halo supportable by Lyman a. 

V. Presentation of a first order infall scenario and the time 
dependent properties of the system it would construct. 

VI. Evaluation of disk properties at the epoch at which 
further material is supportable against infall by Lyman a 
pressure. This includes determination of the two dimensional 
family of disk galaxies which can arise from the three-parameter 
set of perturbations in the Hubble flow. 

VII. Consideration of the fate and subsequent role of a halo 
supported by Lyman a. 

VIII. Presentation of a summary, caveats, promise, and 
directions needing further effort. 



The four frustrations with which this introduction began are 
dealt with as follows: the feedback problem is investigated only 
after the NISM structure problem is eliminated by assumption, the 
star formation problem by derivation, and the infall problem by 
naivete. The feedback result turns out to be significant in that 
it supports the underlying thesis of this paper, namely: The 

eventual properties of spiral galaxies can be controlled by self- 
regulating physical processes occurring during formation. 

In systems like those about to be explored, the supernova 
rate would be quite high. These explosions may have considerable 
impact on the state of the NISM, on the star formation rate, and 
on the interaction of the disk as a whole with infalling or halo 
material. In this initial investigation, however, all such 
effects (but one) are completely omitted. The one effect 
mentioned from time to time is the possibility that supernovae 
would maintain hot matrix (i.e., T ~ 5 x 10^K) interspersed 
through the denser HII phase. This inclusion allows for pressure 
continuity in a system with gaseous filling factor less than 1, 
and in addition provides low opacity channels for photon 
propagation. All further consideration is postponed to a 
subsequent but parallel investigation of disk galaxy formation 
with a supernova dominated NISM. 

It should be clear at the outset that the model presented is 
a prototype rather than a final product. It is the simplest of a 
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class envisioned to incorporate feedback, via physical mechanisms 
of interstellar material, into the description of galaxy 
formation. The general features of the class are summarized in 
Figure 1. 

II# The Self-Gravitating Isothermal Disk 

Assuming that infall is taking place gradually compared to 
the vertical dynamical timescale of the disk, we can 
straightforwardly evaluate the disk structure. For systems of 
interest, mass will remain primarily gaseous during the period of 
model applicability, but a factor f , the gaseous fraction, is 

o 

introduced at the outset to maintain a generality needed mich 
later. It is assumed, however, that the stars when they do form 
have a scale height similar to that of the gas (although a 
detailed model would as usual predict a smaller scale height for 
the stars). 

Conditions in an ionized, isothermal, self-gravitating disk 
with total column density a, gaseous column density a = f a, 

o o 

volume filling factor f Q for the HII phase, and average nuclear 
mass m = (1.4/l.l)m are given by 



p = 2nkT 
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where g max is the gravitational acceleration for newly accreting 
material and p(z) is the horizontally averaged total mass density 
at height z. Values at midplane are given the subscript zero. 
(The hot matrix phase lends continuity to the pressure with no 
contribution to the mass density.) 

The star formation timescale at midplane will be designated 
t m and is given by 
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It is clear that the instantaneous recyclying approximation is not 
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threatened for n < 100 cm , since t is then longer than 3 x 10° 
o ra ° 

years, but note the inverse quadratic dependence on a* 

Let us now explicitly relate the recombination rate to the 

total eraissivity n of the disk, in starlight* This emissivity 

arises primarily in the same massive stars that supply ionizing 

UV. We assumed in the Introduction (following Paper I) that 30% 

of the luminosity of massive stars emerged beyond the Lyman 

limit# The average energy of these photons is <hv> ~ 20 eV and 

each such photon balances one recombination* Thus n is related to 

2 

the recombination rate 0.3 n/20 eV = an , or 


n « (67 eV) an , 


i.e. , the starlight production rate is 67 eV per recombination. 

One might wish to round up to 70 eV or 100 eV as time goes on, to 

include the accumulation of lower mass stars* 

The direct Lyman a emissivity of the gas due to recombination 

of hydrogen to levels with n > 2 is (10.2 eV) ann^ or 
2 

(9.3 eV) an • This may be augmented quite a bit by collisional 
excitation, since the remaining energy of ejected electrons 



(roughly 5 eV on average) is most likely dissipated by excitation 
of Lyman a in a gas with low "metal” abundance. (At least 4 eV of 
the 20 eV available per ionization must however emerge in Balmer 
and higher continua and lines.) The total Lyman a emission per 
recombination (of either hydrogen or helium) can thus be written 
(3/4) 1^ ty, where the factor ij>, the mean number of Lyman a photons 
generated per recombination, is such that the product lies between 
9.3 and about 15 eV. Then 

n(Lya) = In ♦) an 2 . 

The luminosity per unit area, i (total) or l a (Ly a) follows 
from integrating these eraissiviti.es over z. The key ingredient 
from the disk structure is 

/ n 2 (z)dz = n 2 H. 

3 ° 


Recalling that 67 eV » <hv>/f 2 . 


i = 
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Similarly 
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Notice the extreme sensitivity of the surface brightness to a. 
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The mass-to-light ratio is proportional to a • Two further 
things are noteworthy: The uncertain f Q has cancelled out, and 

the Lyman a luminosity is a significant fraction, perhaps 20%, of 
the total. 

The Paper I NISM model thus unambiguously specifies the 
structure, star formation rate, and luminosity per unit area as 
functions of the disk surface density. We are now prepared to 
search for consequences of those properties. 

III. Radiation Pressure in the Disk 
A. Initial Survey 

One searching for radiative feedback is naturally interested 

in properties of the disk when it reaches a column density 

corresponding to the Eddington limit, that is when 

A £/2c =» mg = 27TGom, where A,^ is the Thomson cross section, 
i max r 

_2 

The result, however, (a ~ 0.4 gm cm ) is beyond the range of 
applicability of the model. (The instantaneous recycling 
approximation fails.) A somewhat greater opacity is offered by 
neutral hydrogen (absorbing Lyman continuum) but the limiting 
surface density is still too high to be of interest. 

Another interesting possibility is that a disk emitting 



primarily 0 star radiation might behave like an 0 star itself, 
with a radiatively driven wind* A single isotropic scattering of 
each emerging photon would provide a net outward force per unit 
area (a radiation pressure field distributed through the 
scattering volume like a body force)# 
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Such a pressure would become important to disk stability when 
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comparable to p Q = ttGo f^/2. This occurs when 
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Clearly with f ~ 1/2 and f ~ 1, one would arrive at a familiar 
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looking galactic surface density and scale height at just the 
epoch that the galaxy becomes highly repulsive. 

A detailed mechanism for photon scattering is absent from the 
above picture. However, since a significant fraction of the total 
luminosity is converted to highly trapped Lyman a, whose 
scattering one can hope to model, we now proceed to that specific 



example. 

B. The Optical Depth to Lyman a Scattering 

Assessing the dynamical importance of Lyman a requires 
evaluation of its radiation pressure gradient. In a highly 
scattering system, radiation pressure is essentially isotropic and 
equal to one third of the energy density. Energy density is in 
turn approximately the product of emissivity and average photon 
trapping time. Finally, the trapping time follows from solution 
of the diffusion problem, whose input must be system properties 
including scattering optical depth and line width. 

In order to estimate the optical depth for Lyman a scattering 

we must be somewhat more specific about the distribution of 

material and its ionization fraction. A Strbmgren radius (at 

densities of interest) is roughly a factor of five smaller than 

the scale height H, so that the optical depth of a half plane is 

roughly that of five radial lines of sight through HII regions 

plus any intervening neutral material. An HII region has a 

neutral fraction of order 10~^ and a Lyman continuum optical depth 

of roughly unity. At 10^K the optical depth of Lyman a at line 

center is 10^ times the continuum optical depth, i.e., 10^ per HII 
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region radius, or 5 x 10 for the HII region portion of the line 
of sight in the half plane. This optical depth through 
photoionized regions, although- large, is utterly negligible 
compared to the potential contribution of neutral material between 



HII regions, unless that material comprises less than 10 ^ of the 
total gaseous mass. 

Thus we must attempt to estimate the ratio of interstitial 

neutral gas to gas within ionized regions. (We have already 

assumed in the Introduction that the ionized gas contains most of 

the mass, if not the Lyman a opacity.) A formal estimate of the 

amount of neutral material follows from the requirement that it is 

the location of star formation responsible for keeping the rest 

ionized. We can write the total star formation rate as M /t 

g m 

where M is the total gaseous mass and t the star formation 

g m 

timescale derived in the Section II. Equating this required rate 

with the available rate written M /t , where M is the neutral 

n n’ n 

mass between HII regions and t n is the mean time for neutral 
material to convert to stars, the interstitial mass fraction is 
then 
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where the NISM model tells us t - 4.3 x 10 I years/n_ (cm ). If 

m J o 

gravitational free fall were responsible for star formation, then 
t n ~ t^£ = 4.3 x 10 7 years/[n Q (cm 3)jl/2 Spitzer 1968) and 

M n /Mg ~ [n Q (cm **)/ lOOl^]*^. The implied neutral fraction is 
fairly high at all densities of interest (i.e., 8 to 30 cm as we 
shall see). More likely, however, neutral material is compacted 



by HIT region expansion, stellar wind ram pressure, and supernova 

explosions so that t ft may be roughly equal to the evolution time 

of nearby massive stars. In that case, t^ ~ 4 x 10^ years and 

M n /Mg ~ n o (cm”^)/100I. For the highest densities of 

interest (^30 cm with I ~ 3), the neutral mass between HII 

regions could still be almost 10 percent of the total. In general 

the implied Ly a optical depth ranges from about 10^ when Ly a 
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first becomes important (Section IV) to 2 x 10 . Finally, as an 

extreme lower limit to t , consider that at a density of order 10 

n 

cm , the radius of a region containing 100 is of order 4 pc. 
Its crossing time at a velocity of 10 km s * is 4 x 10^ years, 
making this approximately the minimum time to crush such material 
into a star. With c n ~ 4 x 10^ years, the neutral mass fraction 
is reduced by a factor of 10 from the last estimate above, making 
the highest optical depth of interest 2 x 10^. 

The last estimate, however, brings up an additional 
complication. Although on average each StrSmgren sphere must see 
the formation of 1 massive star during its lifetime, with such a 
short t^, only 1 in 10 would be driving such formation (on a 
fraction of its periphery) at any given time. The distribution of 
neutral material would thus be extremely patchy and it is possible 
that the patches would not get the opportunity to interact with a 
large proportion of the diffusing Lyman a. 

The results of this section are admittedly not very 


satisfying. We are left knowing only that t q , the half plane 

optical depth at line center for scattering Lyman a, probably lies 
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in the range of 10 to 10 with a cross section profile 
corresponding to roughly 10^K. Fortunately this knowledge turns 
out to be sufficient to first order, and we now turn to the 
diffusion problem at high optical depth. 

C. Lyman a Diffusion 

The diffusion of a Lyman a photon out of its parent cloud and 

the disk as a whole will consist of a large number of scatterings 

with short mean free path (involving little time and little 

progress) followed by a small number of scatterings in the wing of 

the line, with larger mean free paths (much more time between 

scatterings and much more progress) (e.g. , Osterbrock, 1962; 

Walmsley and Mathews, 1969; Adams, 1972, 1975). The total 

trapping time is approximately linear in the scale, and 

fortunately the ratio of trapping time to scale is only weakly 

dependent on the line center optical depth t q . Adams (1975) 

calculated trapping times versus for the case of a homogeneous 

plane with photons originating at raidplane. The results (shown in 

his Figure 1) imply typical photon paths of about 15 H (trapping 

times of 15 H/c) Tor T < 10^ and paths of 13 H (t'/IO^)*^ for 

o ^ o 

£ * 

T o 10 • The inclusion of extreme inhomogeneity in both source 

(recombination) and scatterer distributions, velocity structure, 

and the dispersed distribution of the original photon source 


throughout the disk can be expected to lower the average trapping 
time somewhat for our galaxy disk model. The effects are small, 
however, the bulk of progress being made far in the line wing 
where detail is insignificant. Considering Adams 1 results as an 
upper limit and recognizing the extreme insensitivity of the 
result to t q , the trapping time written 

t.iSi 

c 

gives an adequate representation of the situation. The fudge 
factor 6 is expected never to deviate by more than a factor of 3 
from unity. 

We will soon be considering a situation for which Lyman a 
pressure is comparble to or even greater than thermal pressure. 
Krolik (1979, 1983) has suggested that Lyman a pressure in this 
regime will drive an instability tending to fragment the material 
down to a scale such that the line center optical depth is only 
moderate per fragment, releasing photons more readily into 
optically thin crevices. A hot matrix generated by supernovae may 
also provide macroscopic channels for photon escape. At the 
smallest scale, the crevices increase the photon mean free path 
very little if the fragment filling factor remains large. Only if 
channels at each scale lead sequentially to ones of larger scale 
can this raicrocracking lead to an increased escape rate. A rough 



approximation to such a system is one dimensional diffusion along 
a line for which mean free path increases exponentially with 
progress achieved. Taking the smallest and largest mean free 
paths to be H/t q and H/2 respectively, I estimate the trapping 
time to be roughly q £n(T Q /2)H/c where q is the number of 
sequential steps in the same direction between doublings of the 
mean free path* For t q > 10^ and q > 1, t > 13 H/c. Thus crack- 
augmented diffusion seems not to invalidate the previous 
conclusion regarding the trapping time* On the other hand, in 
enhancing escape of radiation close to line center it could 
significantly alter the frequency profile of emerging flux 
(backfilling the very dark line center found by Adams), and hence 
its ability to interact with low velocity material at high z* 

D. Conditions in the Disk When Lyman a Pressure Becomes 
Significant 

The Lyman a pressure at midplane, one third of the energy 
density, is (for a trapping time 15 H6/c) 

/T . 1 * /T n 15H6 156 _ Z 

p(Lya) = j * f Q n o (Lya) • = JT * a ~ 3d 

Since the result is roughly three times the previously estimated 
total single scattering radiation pressure, we can expect it to 
become important at an interesting (i*e*, galaxy-like) column 
density* Other sources of opacity are not likely to alter the 



total radiation pressure significantly. 

The radiation pressure becomes important to the disk 
structure when p(Lya) reaches p Q and the parameters of the disk 
when this occurs are referred to as "critical** in subsequent 
sections of this paper. Characteristic values are given the 
subscript 1. At the critical epoch, the disk parameters are: 
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Expressing the recombination coefficient a in fundamental 
constants, the critical surface density can also be written 
(Rybicki, 1983) 
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where a is the fine structure constant, r the classical radius 
r s o 
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of the electron and N the “large number" e /G m ra ~ 2 x 10 J # 

p e 

The critical epoch at which Lyman a pressure reaches p Q is 
thus characterized by a column density whose magnitude depends 
only on m^, e, and c. It knows nothing of perturbations in 
the universe, galaxies, or even gravity. The scale height of that 
material depends also on TT and G in such a way that it is of order 
100 pc, also with no prior knowledge of the nature of galaxies. 
Both results are sufficiently attractive that it is tempting to 
suppose that Lyman a pressure simply halts all further infall when 
this state is reached. The situation is more complex than that. 



however, and its discussion is postponed to the next two sections. 

E. Three Details 

One may wonder about the simple equation of p(Lya) and p Q for 
determining the critical epoch. In examining Adams 1 (1972) graph 
of the source function versus optical depth, it seems reasonably 
convincing that the radiation pressure gradient Is of order 
p(Lya)/H. In the self-gravitating plane, comparing this gradient 
with Pg is equivalent to comparing p(Lya) and p Q . 

The presence of a small amount of dust in this nominally 
pristine system could seriously deplete Lyman a pressure and 
invalidate the results. This difficulty is present only over a 
narrow range in dust abundance, however, since making the disk 
optically thick to continuum radiation subjects it to the full 
single scattering radiation pressure Z/4 c, effectively reducing 5 
only to 1/3. 

The Lyman a absorption cross section on dust, at solar 
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neighborhood abundances, is roughly 10 cm per hydrogen atom 
(Savage, 1983). Since the Lyman a path is 15H, its traversed 
column density is roughly 15 f a/2 = 7.5 a^, approximately 5 x 
10 atoms per square centimeter when Lyman a first becomes 
significant. Dust absorption is therefore negligible at this 
epoch only if its abundance is less than about 10 of that in the 
solar neighborhood. For a dust abundance greater than 10~* solar, 
the continuum opacity in the optical is important at a * 0.05 gm 



cm , the critical column density found for single scattering. 

The amount of heavy element contamination present in a 

galactic disk when it reaches depends on the length of time 

infall requires to provide that column density and the elemental 

yields of stars with lives shorter than that time. Even with a 

modest elemental enrichment however, it is unclear whether UV 

absorbing dust would have time to form. Suppose for example that 

it forms in the atmospheres of red giant carbon stars. Throughout 

the remainder of this paper, it will be assumed that infall was 

2 7 

rapid (compared, say, to t ral ~ 516 x 10 years) and dust 
insignificant. The converse situation may well merit further 
investigation. 

Notice finally that as a approaches 30 , Lyman a would have 

a high probability of escape via incoherent electron scattering 
into the far wings of the line (Adams, 1983). With these caveats 
in mind, let us now turn to the effects of the Lyman a trapping on 
the structure of the disk. 

F. Disks With Higher Column Density 

For column densities above the critical value, the calculated 
Lyman a pressure in the initial model exceeds the thermal 
pressure. Two competing possibilities emerge. Either Lyman a 
pressure is able to rearrange material such that escape is 
enhanced, leading back naturally to equipartit ion, or it finds no 
rearrangement that will substantially reduce the escape rate, in 
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which case Lyman a assumes responsibility for support of the disk, 

thermal pressure becoming insignificant. In either case, one 

expects a Ray leigh-Tay lor like instability but of a significantly 

modified character (e.g. , Krolik, 1979). The two chief 

differences from the classical situation are that the lighter 

photon fluid is continously resupplied inside the densest portions 

of the denser fluid, and there is considerable interpenetration of 

the two fluids. The cracking described in the last section is 

expected, but the weak gradient in pressure does not lead to 

stable condensations . Lumps continually form, dissipate from 

their own Lyman a emission, and reform, the system behaving as if 

opalescent. If the earlier calculation of crack-augmented 

diffusion is correct, the net escape rate is altered very little 

and a disk supported by Lyman a is the inevitable consequence. 

Rather than do the radiative transfer problem in a radiation- 

supported disk, I shall simply continue to equate p(Lya) to the 

2 

midplane pressure, still necessarily p = ttGo f /2. Other 

O 

surviving relationships are o = 2Hp Q an< * n = fgp/f Q m. In 
calculating the total luminosity, I shall also assume that 

oo 

r 2 2 

J n dz a (4/3) n continues to be a reasonable approximation. 
—00 

By continuing the assumption that the Lyman a trapping time is 
15 6 H/c, while recognizing that a slightly lower value of 6 may 
be appropriate for a radiation dominated plane, the properties of 


the Lyman a supported disk are then 
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The applicability criterion is x > 1/f . The net result of 

o 

Lyman a pressure is to stabilize the scale height so that 

thereafter the density Is simply proportional to q» The rate of 

increase of Z and decrease of t with o are also lessened* The 
a m 

above equations assume that thermal pressure is negligible* The 
actual ratio of thermal to total pressure is approximately x^*, a 
nicety which will continue to be neglected, even though we shall 
find that the validity of the model is severely strained for 
values of x greater than about 10. 



IV. Radiation Pressure In the Halo 

The Lyman a luminosity per unit area of the disk, in 
equilibrium, is Z a - (4c/15S) p Q . If the scale height is 
perturbed from its equilibrium value, is inversely proportional 
to H while p Q is unchanged. Thus Lyman a pressure is able to 
raise the bulk of material in the disk only to a height H^. As 
more material is added, nothing about this conclusion changes. 
Increased radiation pressure will not blow away an excess surface 
density (excess compared to a^). 

A separate question is whether Lyman a pressure can blow away 

material which is not yet a part of the disk. Consider, for 

example, the possibility that Lyman a is. retrapped in a planar 

layer of scale in material with very large velocity 

dispersion. In that case, the already broad profile will quickly 

be redistributed over the effective doppler width of this layer 

and must diffuse to get out. The trapping time in that case will 

be roughly 2 Znx^z^/c (t^ < 10^) where is the effective line 

center optical depth of the layer. The trapped energy per unit 

area (per side) is Zjj/c an< * the pressure p^ 

= l fcnx /(3c) '« 4 p fcnx /45S = f og J,nx„/45<$, where f and a 
ot cl oH g °niax H g 

refer to the values In the disk* The radiation force per gram is 
P H^ P H Z H 2=5 exceeds g max only for Aa < f^a £nr H /A56, 

rather less than' o itself but independent of z^. 


Thus a few 
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percent of a can be undergoing active expulsion at any given 

time. The cumulative effect, however, depends on the time 

required for expulsion. The maximum amount of expellable material 
* 

in time t is 
ra 


a 

e 



f at 


g 

m °max 
456 v 


e 


£nT H 


where v q is the escape velocity. The product at^ - a l t m l 

constant while g = 2 ttGo. Thus 
°max 

a 2TrGa.t t Ant,. 2.3 I £nx 0 km s * 
e 1 ml H H 

a 456 v v 

e e 

Clearly the ratio of expellable mass to total disk mass is small* 

A similar conclusion applies to the stopping of rapidly infalling 
material. 

As we shall see, what the Lyman a pressure can do well is 
support (rather than expel) a rather considerable halo of material 
in the vicinity of the disk. Rapidly infalling material will 
enter, but the much larger amount of material in the vicinity, 
which would enter in time were it not for Lyman a pressure, can be 
held off for a time of order t • Lyman a can also considerably 


*The mass conversion timescale is the length of time the Lyma^n 
alpha can remain bright. 



influence the state of that material, drastically altering its 
eventual fate* 

To illustrate the above point with the simplest possible 
model of Lyman a trapping in the halo, assume merely that once the 
Lyman a has diffused out of the disk (with trapping time 156H/c, 
it is further trapped in the halo* There, owing to the spherical 
geometry, the trapping time is written in terms of R (rather than 
H) and reduced by about a factor of 2, i*e*, it is 156 H R /(2c) 
(Adams 1975), where a new fudge factor 6^ once again allows for 
later improvement. In particular I assume that an annulus of the 
disk, from R to R + dR contributes energy 

15R 6 

dU = -x — - 2ttR£ dR 

ZC CL 


uniformly to the volume interior to R, making its contribution to 
the pressure there 


, 15 6 H 

dp = 4“ F" * 


dR 
a R * 


Similarly (but not quite consistently) the annulus contributes an 
energy 


156 „ r 

2ttR£ dR 

a 


dU 


2c 
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Co Che volume wichin r > R making Che pressure conCrlbuClon aC r 


dp 


15 6 H /R\2 „ dR 

rri? L *ir 


AfCer InCegraCing from R = 0 Co r e , the disk radius, and 
dif ferenCiaCing wich respect to r, the pressure gradient Is 


9p _ 15 , 1 /- r ei3 

9? " " T 6 h T — 

e 


and 

3p _ 15 , l o 1 

Ir ~ ~T Tv 


r > r 


e 


r < r 


Roughly speaking, the energy density is enhanced over the free 

streaming value by the factor 155^/2 at all radii. 

The supportable mass density in the halo is given by 
o 

-3p/3r « pg « pGM(r)/r . For r < the implication is that the 

halo mass within r is directly proportional to the disk mass 

2 

within r. The disk mass is = -nv a, the halo mass can be 

2 

written as = $rcr a, and the halo density is 
2 

P = (dMy/dr)/4irr = go/(2r). The hydrosCaClc condiCion is Chen 


r'.ri-'S " 1 + s) • 3 


or 


3(1 + 3) = 


15 6 H £ a 5 H f g 
211 Go 2 C 6 * 


Thus for 6 f /6 “ 1, 3 ° 0.62; the supportable mass in the near 
H g 

halo (r < r g ) is some 62 percent of the disk. mass. 

Extending the analysis beyond r = r g , we have 


15 . 1 f r ei 3 _ 1 dM 

r'irrir 1 '3 < *i7 

e Att r 


which yields 


m 2 - d + e > 2 ^ + 2 ^ £ f g -i) 

e 

2 2 2 

where Mp = 0 • For r » r e » the system mass increases 

linearly with r. The rotation curve is flat. For example, 

assuming 5 f /6 « 1, g = 0,62, the rotation velocity at r = r 

(neglecting the deviation from a spherical potential) is v rQt = 

1.27 (GM /r )^ 2 while for r » r it is 1.19 (GM /r 
v D e e v D e 

The details of the above analysis depend on the value 

of £ for a radiation-dominated disk. For surface densities 
a 

below Op correct results are obtained by using £ q for that 

regime. At any surface density Lyman a is capable of supporting 

some material for a time of order t . For the radiation-dominated 

m 

disk, the amount of supportable material within r > r e is * 



essentially M^r/r^ 

In discovering that very little material is either expellable 
or stoppable but that a significant amount is supportable, we are 
led to the supposition that infall will continue until such time 
as the amount of nearby material in the infall drops to the 
supportable halo of the then existing disk# Support begins, 
infall stops, and for the most part the material in the disk turns 
into stars in a time of order t m> after which the Lyman a 
brightness drops precipitously and support of the halo material 
vanishes# 

When halo support fails, the material can be expected to 
resume its infall. If as discussed in VIII, however, that 
material has been forced into dense condensations by Lyman a 
pressure during the brief period of support, it will remain as an 
interpenetrating halo rather than join the disk. Assuming that 
this occurs, we can proceed to an infall scenario with a specific 
mechanism for stopping infall at a particular epoch# 

V# A First Order Infall Scenario 

This investigation is concerned with the possibility * that 
individual galaxies suppress their own further growth after 
reaching total masses nuch smaller than the mass of the 
perturbations out ~of "which they are growing# The- ideais-not that— 
galaxy formation is necessarily inefficient, only that by * 



suppressing further growth of dominant perturbations, slower 
systems are given time to develop. 

A nearly spherically symmetric (and therefore improbable) 
density enhancement growing in a normal matter dominated universe 
can be characterized by three parameters; the minimum density p ^ 
reached by the central peak before it begins recollapse, the total 
bound mass Mp of the perturbation, and the external torque or 
velocity shear to which it is subjected. The minimum density 
parameter controls the time of initiation of "galaxy formation," 
the mass parameter governs the flow rate of material into the 
"galaxy," and the shear parameter controls the spatial extent of 
the forming disk. 

Concentrating on material very close to the peak of the 
perturbation (the material the galaxy will eventually admit), 
conditions are particularly simple. Galaxy formation commences at 

1 /o 

time t = [3n/8Gp mln ] ' ABB, with core formation. At half that 

time the mass was at minimum density, p . , and maximum 

min 

extension. Taking the shear as linear over the central portion of 
the perturbation, the angular momentum distribution (after 
azimuthal sharing during collapse) will be that of a uniformly 
rotating sphere. Returned material will thus be broadly 
distributed over the disk with a maximum angular momentum per unit 

mass j proportional to the square of the maximum distention 

2 /3 

r ^ • Hence j « M ' where M is the returned mass at any epoch 
max max J r 



and equals (4/3) P||lln . 

By following the perturbation growth at early times, in both 

velocity and density, and assuming that near the peak the density 

drops quadratically in radius, the return time for the shell 

containing mass M can be calculated* The result for short 

3/2 

times At = t - t after initiation is M c M D (2At/t ) . Thus 

C r c 


j « M 2 ^ 3 « r 2 (M) « At. 

J raax max 

The disk formed by time At has mass M and maximum angular 

1 12 2 

momentum per unit mass i ^ = v . r « [GMr 1 <= r (M) where 

r J max rot e v e J max 

v^ ot is the rotation velocity and r^ the disk edge radius. Since 

M * r , this relation for 1 „ necessarily implies that r ./r 

max max g rnax 

is constant during the infall. In addition, r /r « 
a * e max 

^max^ r max so t ^ at recontract i° n ratio is proportional to 

the square of the shear. With this result it is convenient to 
represent the constant shear parameter by e defined such that 


r 

max 

Since r (M) * (At)*^ while M « (At)^^, the disk radius r and 
max e 

1/2 

average surface density a are each proportional to (At) so long 
as 2At « t c (equivalent to M « Mp, the case of interest for a 
self-limiting galaxy). 

The characteristic parameters of the unimpededly developing 


disk as functions of its mass are then: 



r = (3M/4np ) 1/3 = (2MGt 2 /ir 2 ) 1/3 

max x min c 
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t = 2 tt t f1 = 2 3/2 e 3 t 
rot flow c 


o = o/ (2At). 


Notice that the flow timescale (or rotation period) is a 

3 

constant, proportional to e t^ for the forming galaxy. From the 
point of view of the galaxy itself, this parameter combination is 
one of two which govern the eventual structure. The galaxy scale 
at any value of a follows directly from it: 


r 

e 


ir G o t 


2 

flow # 


The other significant parameter is the length of time required to 



reach some particular value of o, and in the next section At^ is 

introduced, the time to reach and the Lyman a domination of the 

disk. The three parameters of the infall reduce to these two 

which govern the nature of the galaxy formed. 

It is of some interest to estimate the value of e that would 

be appropriate to the Milky Way Galaxy, checking on the internal 

consistency of these results. Assuming a mass and formation time 

11 9 

respectively of M ~ 10 and ^10 years implies * max ~~ 45 

kpc. It appears that the mass in our galaxy never extended 
farther than the Magellanic Clouds. In addition, for an edge 

radius now of 10 to 15 kpc, e is roughly 1/2. That translates to 
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a rotation period of 3 x 10 years or flow time of 5 x 10 

_2 

years. The corresponding value of a is roughly 0.07 gm cm , 

300 pc or 5 a^. 

It was implicitly assumed that material flowing into a galaxy 
joins the disk quiescently. This point of view is not altered 
drastically by a cooling time estimate for shock heated infall 
material. If a halo containing mass AM < M within radius r is 
heated to T ^ GMm/(5kr) and has a cooling time t C0Q ± - 3kT/Ln 
where L is the bremsstrahlung cooling coefficient and n is the 
average halo density, then it can be shown that 


t cool 
C f low 


M r 
AM R 

coo 


1 



where 


R cool 


°- 95 “£, ( w -> 


Vl/2 


Gm m o 
P e 


35 kpc. 


where a- is the fine structure constant and r the classical 
f s o 

radius of the electron* Thus near halos with AM ^ M will cool in 

less than a flow time by bremss trahlung alone for galaxies less 

than 35 kpc in radius. Note that the Lyman-a-dominated-disk scale 

height can be written » 24a,. (m /ra )*^ 2 R ' = 4 x 10” 3 R 

° 1 fs e p cool cool* 

that the cooling radius is independent of the mass within it, and 
that its value is disquietingly similar to the value of r max found 
for the Milky Way. 

One might also be concerned about the amount of star 
formation which occurs in the infalling material before it reaches 
the plane. Such star formation will tend to be suppressed by 
ionizing UV leaking from the plane. But even if not suppressed, 
it is small. The infall is not as dense as the plane, and the 
infall time is less than the age of the plane. If they both have 
the same density dependence for their star formation rates, 
infalling material will always be much less stellar than the 
material in the disk. 

A bit more detail about the early behavior of the growing 
perturbation may be of some interest and is included in the 
Appendix. The next section considers the properties of disks 



whose orderly infall is interrupted by the onset of Lyman a 
support in the halo* 


VI. Disk Formation Regulated by Lyman a 

At the epoch when a reaches c^, taken as fiducial, the system 

2 2 3 2 3 4 

radius is r. = ttGo. t-_ the mass M. = ira.r, = tt G a. t-. , 

1 1 flow 111 1 flow 

and the incremental time 


. . 3 2 3 4 

Aj _ c r M ^2/3 c °1 flow^2/3 

* r i — b; > • 


Expressed in terras of an( * Atj conditions at other times are 

simple functions of 


x = o/a 1 = (At/At^ 1 ^ 2 . 


The gaseous fraction, f , as a function of a must be 

o 

evaluated before the mass of the supportable halo can be 
calculated* The UV photon production per unit area is 
(4/3)Hf o cm o 2 so that the rate of increase of stellar column 
density (for o > o^/ f ) is 


o A = — ~ ( 1 . A m ) 4 - Hf an 

* m 6 p 3 o 
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— T— f Han. x 
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Using the relations o = a/( 2At) and 2At = 


have a = o^/(2At^x) and 


da * °* 2 2At l 3 e 4 At l 3 ( ° ~ ° 
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Substituting 


1 r 4 At l^ 3 
T hr t — x , 


J -i(T rf 


the required integration can be written 


x, e -J(x) J^x) 1/3 J 

— = 1 - f = r-pr J J e dJ. 

x 8 [JCx )] 1/3 o 


When integrated numerically and fitted, the result is 


g 4 + 3J + J 2 
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This result together with the definition of J provide the 
dependence of f on o during the unimpeded infall period* (It was 
assumed, however, that f^x > 1.) 

In §IV we found that the ratio of supportable near halo to 
disk mass, 3, follows from 3(1 + 3) = f^/6. We now neec * on ly 

to examine the infall to discover when the amount of nearby 
infalling material becomes supportable. This can be done in 
sophisticated ways, but in the following approach the essential 
elements are clearer and the dependence on infall details is 
weaker. 


Nearby infalling material would enter the galactic disk over 
a time period comparable to ® r e / v ro t* Hence the amount of 

material trust be of order Ao *=* a t flow = 0 t flow /(2At) and 
becomes supportable when Ao » 3o. Thus support commences when At 


reaches For 3 ~ 0.6, unimpeded infall lasts only for a 

time of order 0.8 tf^ ow ~ t rot /8, roughly 4 x 10^ years for the 
Milky Way Galaxy. 


The disk surface density at infall shutoff follows 
from o = Oj(At/At^)*^, 


x 


f C f low\ 1/2 
l 23At x J . 


To find the resulting family of possible galactic structures ,_it_ 

is most convenient to calculate the family loci for constant J. J 



leads immediately to f (J) and 3(f )• 

O O 

Mapping can proceed using x as the implicit variable, 

“l ' I C ml J/x3 

2 9 

At - x Atj - T t ml J/x 
^ f low = 2BAt = T 6t ml J/x< 

When such mapping is done, the parameter pairs (t^ Qw , At^) for 
J > 1.5 are found to duplicate (at larger values of x) the values 
for J < 1.5. This makes J > 1.5 (f < 0.37) inaccessible and over 

O 

most of the range available to Atj) the simpler case f^ » 

1 suffices. 

Before performing an actual mapping, it is useful to discuss 

2 

the rotation velocity. We found previously that r = 7r( ' at fi ow * It 
is unclear in this simple model whether a should include the 
supported halo but I will assume that it does. Thus 
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The third-to-last form for v roC is particularly illuminating 
in that it reveals the complete set of assumed parameters leading 
to the final rotation velocity of the galaxy. Lumping all of the 
uncertainties in fjf 2 ^3 ^ nto * (see §1), the only other parameter 
of significance (given J, which determines fg) is 6^. 

In addition, f J is monotonic in J for the applicable range 

O 

J < 1.5 so that the maximum value of f J is 0,558; the maximum 

O 

possible rotation velocity is 129 I 6^ km s Clearly our 
mapping of galaxies in the At^ plane will produce familiar 

results only if 16 >2. A slight weighting of the ZAIMF toward 

higher masses, compared to the solar neighborhood IMF, is 



45 


sufficient and perhaps expected. Once again the model produces 
galaxy-like numbers without prior knowledge. For purposes of 
illustration, 6^=1, 6 - 1, I 3 3 will be assumed in the 
numerical example below. The largest rotation velocity will be 
388 km s 

Figure 2a shows the properties of disk galaxies for which 
Lyman a pressure stops infall, as functions of the two composite 
parameters t^^ Qw and At^. Recall that t^^ Qw is 1 / 2tt of the edge 
rotation period (constant during the formation of a disk) and At^ 
is the time required for a to reach The infall duration At is 

tfio W /(20) where 3 varies between 0.29 and 0.61 over the diagram. 
Hence the timescale for accumulation of the "galaxies” shown 


is At ~ t 


flow* 


Figure 2b presents the same family of galaxies as functions 

of parameters characterizing the perturbation as a whole. 

g 

Assuming that galaxy formation took place at t c ~ 10 years, the 
figure presents 
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where = tt(ttGo^x t^^ Qw ) o^x is the disk mass at the time support 

Q 

begins. For t = 10 years, Mp’ = Mp is approximately the total 

2 

perturbation mass, e' = e is the shear parameter, and £ = r e^ r max 

is the net shrinkage. Values of tf^ ow ** At which are small 

Q 

compared to 10 years imply moderate to small e and galaxy masses 

3 / 2 

small compared to M p . The fraction M/Mp is in fact (2At/t c > ~ 

^ and is shown along the right vertical axis 

in Fig. 2b. 

In both representations Figure 2 locates galactic disks on a 

bent triangular sheet showing loci of constant r e , o, and v rot « 

Along the rolled edge the forward (accessible) section is bounded 

by a line labelled by* the maximum rotation velocity 388 km s * (J 

= 1.5). Beyond this edge, no Lyman a stabilized systems exist. 

The other two edges of the sheet are chosen for convenience. The 

-2 

edge labelled 0*015 gra cm is the limit of disks which can be 
supported primarily by Lyman a. Similar models can be constructed 
beyond this edge using the Lyman a emissivity of a thermally 
supported disk. The strong dependence of Z a on a in this regime 
leads to a more gradual variation of a than that shown within the 
triangle, but the basic scheme of the model is preserved. (The 
precise boundary of applicability is shown by the dotted line 
along which f^a = a^.) Dashed lines in Figure 2b show extensions 
of the r = 0.25 kpc and r « 1 kpc loci beyond this edge, as well 
as the x = 0.5 (a = 0.007 gm cm ) trajectory. 



The third boundary is labelled r = 0.25 kpc. Solutions 
beyond this boundary are also possible in principle but the 
calculations involve spherical symmetry since r and H are 
comparable. Preliminary calculations indicate that the 
characteristic scale R c remains about Hp depending on the filling 
factor. The clear distinction between disk and halo disappears 
but there is a change in p(r) dependence from inner to outer 
regions, at about R c * Apart from the inner few hundred parsecs 
such systems have M r but no characteristic mass. They seem to 
be dwarf spheroids although I use the name with some caution since 
they potentially extend to very large radii. During formation the 
inner region would seem to resemble the huge HII complexes seen 
within other spiral galaxies. 

There are two other applicability limits, shown as horizontal 
dashed lines in both diagrams. The upper line marks e f = 0.5 
(tfi ow = 5.6 x 10’ years) and reflects the fact that the maximum 
allowed velocity shear in a gravitationally bound system results 
in net contraction by a factor of four once angular momentum is 
shared azirauthally. Thus the universe is unable to supply 
perturbations with £ > 0.5. Notice that this upper limit just 
barely admits the Milky Way. 

The lower line at = 3 x 10^ years (e* = 0.19) marks the 

failure of the instantaneous recycling approximation to 
accommodate to changes on timescales of order t^ 


or At. A 



related applicability limit is that the star formation timescale 
must exceed 3 x 10^ years. For I = 3 as in this example, t ffl = 

1.3 x 10 9 years/n Q requiring n Q < 400 cm -3 , o/oj < 50, a < 0.7 gm 
cm • Recall that incoherent electron scattering also becomes 
important in this regime* More conservatively, the model's 
details are in doubt over much of the corner shovm with o > 0.23 
gra cm * 

The total range of parameters for which this model (or its 
cousins at low surface density or small scale) could be applicable 
is, on Figure 2b, to the left of v r « 388 km s * between e « 

0.19 and 0.5, excluding the high surface density corner at lower 
right. It thus contains only systems which are rather severely 
sheared, but a very large range in perturbation mass. It also 
contains a very broad range of resulting galaxy parameters. 

Characteristic scales range from roughly 200 pc to 20 kpc. 

—2 

Surface densities range from roughly 0.01 to 0.3 gm cm • 
Furthermore the fraction of the perturbation mass which enters the 
galaxy ranges from a few percent down to 0.0007. 

Despite the relatively high specificity of the model thus 
far, most of the conclusions one might be tempted to draw from 
these results fall in the regime of wild speculation. The problem 
is that when a system with relatively low e forms a "galaxy," only 
a small fraction of the total mass is consumed before infall to 
that object stops, or rather pauses, for a time of order t^. 


The 
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eventual structure of the larger system depends critically on what 
else is going on, in addition to the effects of the pulse of 
Lyman a on the surrounding material. 

Consider, however, the high shear end of the possible range 

of conditions. For e « 0.5, radii range from roughly 4 to 20 kpc, 

-2 

surface densities from 0.015 to 0.08 gm cm , rotation velocities 
from 80 to almost 400 km s~*. All three of those variables alter 
in direct proportion at constant e. Roughly 3% of the 
perturbation enters the disk, a similar amount is held in the near 
halo, and the rest is supported (at least temporarily) in a 
massive halo reaching roughly 10 disk radii. Were it not for 
their needing other nearby mass perturbations in order to achieve 
the high shear, these systems would be good candidates for 
isolated spiral galaxies with massive halos. Owing to the shear 
requirements, however, they are more likely to occur near other 
substantial mass concentrations with contraction times at least as 
short. Such systems, in general, are more likely among the later 
forming objects (larger t c ) for several reasons. By forming on 
cluster boundaries where they began as positive density 
perturbations in comparatively lower density regions, their 
collapse delay subjects them to maximum shear at maximum extent. 

In addition because they form relatively slowly (extended At) they 
compete badly for dominance in regions containing significant 
amounts of low shear material, that is within clusters. 



One other aspect of high shear systems is that the amount of 
halo material with which they nust cope depends more sensitively 
on £ than does any other feature. If such material eventually 
generates the bulge component, the globular cluster host, the dark 
halo, small satellite galaxies, or the outer HI distribution which 
is used to measure the rotation curve beyond the optical disk, 
such manifestations may because of this sensitivity appear rather 
uncorrelated with the gross properties of the disk component. 

VII. Halo Evolution Considerations 

Assuming that halo material is more or less homogeneous when 
it first stabilizes, its properties are defined by the radiation 
pressure available to support it. For r > r e , the density of that 
material is 

* S u i 
H a 

27T cr GM(r) 


where is the total Lyman a luminosity of the disk. Thus 


105 ui 

M(r) - ( — - G ” ° ) 1/2 r. 


The ionization state of this material depends on the amount of UV 

flux leaking out of the plane. The leakage fraction f should be 

roughly 1/3 r (0)/H where r (0) is the typical Stroragren radius at 
s s 

midplane. Assuming the typical 0 star ionizing luminosity to be 



1.5 x 10* 9 photons per second (e.g. , Abbott, 1982) the leakage 
fraction for a Lyman a supported disk is then 
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and the ionizing flux in the halo at radius r 
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With these results, the total and neutral densities can be shown 
to be 
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where A c is the average ionization cross section. Since both of 

2 

these densities are inversely proportional to r , the ionization 

fraction is constant. The residual optical depth beyond r to 

00 

continuum radiation is T c ( r ) " / n°(r)A c dr or 
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For moderate f y> the homogeneous halo is optically thin and fully 



ionized. Its residual optical depth to Lyman a, however, remains 


•i 5 

significant to at least r ~ 10 ~ 10 pc. 
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The halo pressure is 56^£ a y(4;ir c) and the characteristic 


1 / 2 

velocity v = (p/p) is both independent of r and essentially 
H 

equal to v CQt . 


The halo material is thus capable of rearrangement 
with very high velocity and the characteristic dynamical time is 


r v _ = t c , . at least in the near halo. Consistently, the 
e rot flow’ 

3 1/2 

infall time when support fails is t^ n ~ (r /GM(r)) ^ (r/r^) 

tflow* One ot ^ er timescale of interest is the recombination time 


'rec an H,r 
1 e 


1 ^ 1.6 x 10 5 yrs r^ 

an n H,r 

o o 1 e 


Typically t> rec is small compared to t^ Qw at r ^ r g but becomes 
appreciable at large radii. 

The nature of the Rayleigh-Taylor-like instability between 
radiation and the gaseous fluid differs between disk and halo. In 
the halo the radiation source is predominately external, the 
Lyman a arising from the disk. In addition, the severity of the 
overpressure (radiation versus thermal) is much greater, 
essentially (v /10 km s ' versus only x in the disk. These 
differences lead to much greater potential for compression of halo 
material into dense clumps. I have not, however, been able to 
discern a reliable criterion for judging the crossover point 
between the opalescent behavior imagined for the disk and the 
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wholesale clump formation which I propose for the halo. It is a 
question deserving careful study well beyond the scope of this 
initial investigation (cf. Krolik, 1979). 

One can inquire, however, about what sort of clumps might be 

expected to form in the halo under the effects of the extreme 

Lyman a pressure. Such clumps might stabilize when their thermal 

pressure equals the Lyman a pressure; but, that supposition 

ignores the penetration of Lyman a. The latter invades clumps 

easily, and they experience only its relatively mild gradient 

(which supported them before compression). Lyman a pressure is 

not particularly isotropic, however, at the relatively low optical 

depth of the halo. Thus stable clumps will be ones with a 

moderate Lyman a optical depth of their own, experiencing a net 

compressional force per unit area by differential scattering of 

the outward flowing Lyman a component. The outward Lyman a flux 
2 

must be ). It carries momentum equal to that divided by 

c. The fraction reflected by a clump of optical depth unity at 
line center will be roughly 1/x* where x* is the half width (in 
doppler units) of the line, as discussed by Adams (1975). The 
line width is not particularly sensitive to t q so I will set the 
fraction at 1/8, a typical value. We then have the dual 
conditions 


Ap 



_1 

c 




and for clumps of scale R 


V° R - 1 

together with the internal pressure balance Ap ~ 2nkT (with 

4 2 

T ~ 1CT K) and the ionization balance FA n° = an . A is the 

c a 

Lyman a scattering cross section at line center, approximately 10^ 

—2 

A c * These conditions lead to criteria with n <= r , so that the 

compression factor is independent of r# The clump scale increases 

2 4 

as r and the clump mass as r • 

One might hope that these clumps eventually lead to self- 

gravitating. entities and that the near halo produces stars and the 

galactic bulge, the intermediate halo makes globular clusters, and 

the far halo makes dwarf galactic companions# The numerical 

results, however, are not so clear cut. The incremental 

pressure Ap equals only 6/40 of the full Lyman a pressure of the 

a 

halo, or 



where another fudge factor r 6^ is introduced to follow the 
effects of the uncertainty in clump internal pressure# The clump 
density is 
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and the compression factor 


r 5 x f 

ec ( 4 *) 


1/2 


140H 


1 


Only relatively large high surface density systems (r e x > 140 H^) 
seem vulnerable to any compression at all. The neutral fraction 
in the clumps is rather small 
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as are the clump radii 
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Finally the masses of the individual clumps, still for r > r g 
only, are 
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The result for the clurap mass is extremely sensitive to a number 

of factors, making even the order of magnitude suspect* For a 

galaxy like the Milky Way, however, with x * 5, clump densities of 

order 1 cm"^ are suggested at r = r g . The other parameters at r = 

_2 

r g are a flux leakage fraction f y ~ 10 (for f Q ~ 1/2), a clurap 
radius ^ 4 x 10 pc, and a clurap mass ~ 2 x 10 J M € ~ 7 M^. 

Masses would be substantially larger (at r = r^) for lower surface 
density systems. They are also sensitive to the flux leakage 
fraction which depends most sensitively on the formation 
distribution of high latitude 0 stars. Overall the possibilities 
run from substellar to globular cluster masses at least# But if 
the analysis contains even a shred of resemblance to reality there 
should be strong systematic dependencies of clurap mass on system 
column density and, within a system, on r/r e . 

This initial clumping of halo material is not the last word 
in its evolution. As the inner halo forms clumps, it is possible 
for it to become optically thick in Lyman continuum, shutting off 
ionization and heating in the outer halo# The outer clumps can 
then cool rather rapidly and recombine. They may . be induced to 
fragment as their Lyman a optical depth rises sharply. A similar 
effect may apply during the decrease in brightness of the plane as 
star formation lowers the gaseous mass needing to be kept ionized 
in the disk. Once both ionization and Lyman a are effectively * 



shut off, and barring any consideration of effects that supernovae 
would have on the halo, there must be a competition between 
cooling/recorabination on one hand and pressure expansion/inf all on 
the other* If clumps manage to become self gravitating, they may 
individually survive* If not, collisions between them can lead to 
transient high pressures and eventual stability of coalesced 
groups. 

Finally, although one of the ground rules of this 
investigation has been to omit consideration of supernova effects 
as much as possible (in part because those effects promise to be 
as significant as those attributed to Lyman a and deserve a 
separate investigation before a merged mpdel is attempted), the 
relatively low density in the halo makes It a very likely location 
for supernovae to dominate the outcome* Waves of supernova-driven 
compression, followed by cooling and condensation could well 
shatter the halo* 

The above discussion does not leave one confident in any 
particular outcome of halo evolution* But if we cannot clarify 
the answer, perhaps we can at least clarify the problem* What one 
hopes for from halo evolution is that it will generate both the 
observed stellar halo components as well as a largely unobservable 
dark halo of substellar objects* But as we found in the previous 
section, one hopes for this outcome only in the highly sheared 
systems leading to halos some thirty or so times more massive than 



the disk. For less sheared systems, typically with higher surface 
density and smaller scale, the residual mass in the perturbation 
ranges up to a thousand times the mass accepted by the galaxy. In 
that case, one supposes there is some natural cutoff in the amount 
belonging to the "central object" and that the rest is simply 
given time to join other perturbations. This may lead to a 
cluster of galaxies only the central member of which followed the 
formation scenario herein prescribed. 

In regard to this latter possibility, it may be important to 
realize that there is a fourth parameter of galactic infall, 
discussed recently by Silk, Szalay and Zel'dovich (1983). That 
parameter has to do with the non-centrality of the perturbation 
flow. It is a defocussing effect that destines different regions 
in a hierarchical or many-overlapping-wavelength system of 
perturbations to aim for different centers. Such a parameter is 
similar to e in that in the absence of dissipation it leads to a 
system of finite size. If all major perturbations contain similar 
values of this defocussing parameter, and if in general it leads 
to only a modest recontraction, then the low e high surface 
density systems may never have an opportunity to form. One would 
then get disk galaxies only when the shear parameter was large 
enough to dominate over defocussing. All less sheared cases would 
lead to spheroidal systems dominated by dissipation, somewhat like 
Larson’s (1969) models. In short, it is clearly important to 



study halo evolution for the systems with sufficiently high z that 
the galaxy disks look familiar* It is not so clear whether the 
lower e models have any counterpart in reality. If they do, they 
correspond either to a galactic cluster with an unusually high 
surface density disk galaxy (at least initially of small radius) 
at the center, or possibly on occasion to a very massive 
spheroidal galaxy with a high surface density disk buried deep in 
its core. 

VIII. Summary 

The model prexented thus far is largely an impressionisstic 
sketch of how galaxy formation might take place. The individual 
brush strokes in the previous sections may bear little resemblance 
to detailed behaviors they represent. Having put them to canvas, 
however, we now step back for three views in broader sweep. The 
first of these perspectives is a summary of the components of the 
picture; the second is a portrait of the formation of a particular 
system, the Milky Way; and the third is a landscape of the 
possibilities represented in Figure 2. 

The most specific aspect of the model is its assumed 
understanding of the state and self limited star formation rate of 
the nascent ISM. Combining this understanding with self 
gravitation, many properties of the forming disk can be calculated 
as simple functions of the column density o. These include p, H, 



In order to find the last pair. It 


n, n, n , 1 , A , and t or p*. 
is necessary to make an assumption about the high mass weighting 
of the zero age initial mass function (ZAIMF) for stars and the 
parameter I was introduced to express the weighting relative to 
the Paper I value for the solar neighborhood. (The Paper I 
weighting, based on Abbott 1982, was already somewhat heavily 
weighted toward high mass stars as compared to some earlier 
estimates. See Garmany, Conti and Chiosi 1982, for a comparison 
between estimates.) 

Sources of radiation trapping in the disk were explored and 

it was found that Lyman a pressure in the disk becomes comparable 

to thermal pressure as galaxy-like column densities are approached 

(envisioning a monotonically increasing due to infall). The 

competing possibilities that at higher column densities Lyman a 

will either enhance its own escape or take over pressure dominance 

were investigated and the latter held sway. The structure of the 

disk was then recalculated assuming Lyman a support. The chief 

result was that the scale height stabilized at H = 138(6fg/f Q )pc, 

2 2 

where 138 pc = 23.3 (m /m ) a (e /Gra m )r . The crossover 

e p ts p c o 

column density between thermal and Lyman a support is at of 6 = 

6 

0.0146 gm cm -2 = 70 pc“ 2 = 0.28 (kT/I H ) 3/2 m/A T » 0.006 m /Aj. 

—3 

The characteristic midplane density at crossover is 8 cm and the 
mass conversion timescale about 5 x 10^ years. 

* 

A simple model was explored for growth of a perturbation in 



the early universe. In particular, relevant results were 

presented for the infall of the core of such a perturbation. 

Initially, three parameters describing an arbitrary but 

spherically symmetric (except for shear) perturbation were 

introduced, the perturbation mass M^, the initiation time of disk 

formation, t c , and the degree of shear, e. It was shown that 

these lead to two interesting formation constants for the 

subsequent galaxy. The net shrinkage of each mass shell from its 

maximum extension r to subsequent disk extent r is r /r 

max e e max 

2 

= e . The upper limit to possible values is e = 1/2. In 

addition, t flQw = r e v rot - (r e 3 /GM) 1/2 = (2 1/2 /tt)£ 3 t c . Finally, 

1/2 

as infall continues, r g and o are both proportional to (At) and 
M “ M d (2At/t J 3 ^ 2 for 2At « t . 

To this point the overall model provided a time dependent 
picture of a galaxy disk as it accretes matter and assumes an 
evolving structure consistent with its current surface density. 

The next task was to learn when the galaxy was capable of 
rejecting further infall. The mechanism explored was again 
Lyman a trapping, this time in the halo. It was found that when 
the near halo (r < r g ) contained a total mass comparable to what 
was then in the disk, the Lyman a pressure gradient could support 
the halo mass against gravitational attraction. More 
particularly, when the near halo contained A a = So with 0 ^ 0.6 
then Ao was supportable. This leads Immediately to a * 



stabilization time At = t c1 /2 0 (from a = o/2At <• Aa/t,. 

flow flow 

= Thus the duration of unimpeded infall is comparable 

to = r e ^ v ro t* A H- P r °pe r tT es (essentially o, r g and 

combinations thereof) follow immediately. Summarized, they are 
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(This formula gives the disk radius at At* The results given 
earlier included an additional factor of 1 + 8 expected for disk 
growth near the edge where the halo mass was not vertically 
supported* ) The disk mass is 


M 


Mp (2*/g>« - Mp ('L/) 3/2 . 


These results are appropriate within the triangular sheets of 
Figure 2 (or with 3 ** x = o/o^ they extend beyond the low surface 
density edge)* Conceptually similar results can also be found 
beyond the low radius edge, but these systems are spheroidal with 
constant characteristic radius* In principle the above formulae 
should be the most significant consequences of this endeavor, but 
only after the fate of the halo and universe f s choice for e are 


understood 



The "rolled edge" limit of applicability of the above 
equations (in Figure 2) corresponds to systems which collect so 
much material that they turn to stars before stabilizing. Beyond 
this edge (t m < approximately), Lyman a pressure never 

stabilizes the infall and the model as it stands would predict 
that the entire perturbation mass would gradually fall onto an 
essentially stellar disk. Such largely stellar systems are very 
likely candidates for supernova control of disk limitation that 
will be explored in a subsequent paper. If such a model works, 
the resulting galaxies should still be systematically larger and 
higher in surface density and rotation velocity than Lyman-a- 
stabilized disks. 

The rolled edge sets a limit to the rotation velocity of 
Lyman cx-stabilized disks equal to 


rot 


(r V) 

9 6 h JL <dr <°-° 07 hW c 


« 130 I 6^ km s 


It is only in setting the location of this limit that the 
parameter I enters the model results. If 6^ is indeed 1, and if 
ignoring the flattening of the gravitational potential has little 
effect, the value of this limit argues for I ^ 2 to 3 in order to 
reach the observed rotation velocities of galaxies. 



The mass distribution in the Lyman a supported halo is such 
that the rotation curve of the system is flat beyond the edge of 
the disk. In order for this situation to be maintained after 
Lyman a support ceases, the halo mass must clump into long raean^ 
free-path "bits" before infall resumes. The time to do this is 
comparable to the time the disk, at much higher density, takes to 
form itself into stars* Presumably the condensation must be 
driven by energy percolating out of the disk. 

Lyman a pressure in the halo is high, very much higher than 
thermal pressure, but the gradient is weak. The estimated 
clumping factor was not very great and the clump mass calculation 
was subject to enormous uncertainties. A proper calculation of 
halo evolution is the most important next step in developing this 
perspective. It would seem that potential exists for radical 
shattering of the halo, perhaps even into very small pieces. But 
all detail is lacking. 

One point not made earlier involves predictability of 
elemental enrichments in this model. The bulk of disk material 
forms into stars in a time of roughly t m ® 5 x 10^ i/x years. But 
during that time, a fraction * 0.06 I went into 0 stars. Of 
that, a fraction fj » 1/2 was burned to helium or beyond. Given 
that such stars leave a negligible amount of this processed 
material in their stellar remnants after explosion, the enrichment 
level toward the end of this epoch should be 3l % of the initial 



hydrogen mass. Quite possibly the IMF evolved rapidly with this 

enrichment, with an average value of I of order 2 to 3 (from the 

limit on v ). Still one would expect the low mass stars in this 

first generation to be enriched in the sense that some three to 

nine percent of their anticipated hydrogen mass will have been 

converted to helium or heavier elements. The appropriate 

elemental yields will be those of stars whose lives are shorter 

than t • (Recall, however, that when the dust cross sections per 
ra . 

hydrogen atom reaches 0.01/x of the solar neighborhood value, 

Lyman a become subject to absorption. Can Lyman-cx-absorbing dust 
form in 10^ years?) 

Among the specific examples calculated in making Figure 2, 

one which perhaps represents the Milky Way has Mp* = 1.7 x 

12 9 

10 and e f = 0.46. Assuming, as there, that t c = 10 years 

and 1=3, the Galaxy reaches a surface density at At^ = 2.8 x 

6 —2 

10 years and its final surface density 4a^ = 0.06 gm cm - 280 

—2 7 

pc at At = 4.5 x 10 years. This point is reached when f = 
0.7, 0 = 0.47, J = 0.5. Thus the system is still primarily 
gaseous, with a near halo mass equal to about half the disk mass. 

The flow time is 4.3 x 10^ years, the edge rotation period 2.7 x 

8 ““1 
10 years, and the rotation velocity 242 km s • Its nominal disk 

radius is 7.2 kpc but flow of near halo material to the edge 

causes growth to 10.6 kpc. The original disk mass is 2.7% of 

Mp*. Adding the near halo raises this by (t + 8) to 8.5% or 1.4 


x 10“ M^. The reraaiaing material of the perturbation is 
initially supported by Lyman a pressure, arranged with M a r, - * 
providing a flat rotation curve to 100 kpc. Support is 
maintained, however, only for a period of about t m = 4 x 10^ years 
while the disk finishes most of its star formation. Thereafter 
halo material at a distance r from the Galaxy takes about 4 x 10^ 
(r/10 kpc) years to fall in. If it has fragmented into small 
stable bits before reaching the disk, it can pass through without 
dissipation providing a long-term stable halo. The halo will be 
dark today if the bits are predominantly substellar in size. 
Globular clusters and a galactic bulge may derive from star 
formation during the halo fragmentation era. 

The disk itself, when infall is first stopped, has a surface 

—2 —2 
density a = 0.06 gm cm , a gaseous surface density 0.04 gm cm 

22 —2 —3 

(or = 2 x 10 cm ^), a central density n Q - 32 cm , a 

-10 -3 

midplane Lyman a pressure of 2 x 10 dyn cm , thermal pressure 

of 0.9 x 10 dyn cm , a scale height of about 200 pc, a mass 

conversion timescale of 4 x 10^ years, and fc^/c = 5 x 10“ dyn 

cm » Its Lyman a luminosity is 2 x 10^ erg s , its total 

46 -1 

luminosity about 10 erg s , a truly dramatic object! During 

o 

the star formation epoch it burns 2 x 10 of hydrogen, about 4% 

of the disk mass. Much later in the evolution of this system, 

2 3 

when its luminosity has dropped by a factor of 10 to 10 and its 
gas content by a factor of 10 to 10^, its inhabitants will observe 



an apparent equipartition between the starlight energy density and 
the gas pressure. Both are feeble remnants of an earlier time 
when their comparability was important. 

Before leaving the Milky Way, notice that the original star 
formation epoch in this model lasts roughly one sixth of a 
rotation period. The Galaxy could not yet have known about the 
potential for spiral structure. In addition, the galactic plane 
would probably be somewhat poorly defined during this time. 

Locally it is the meeting place of equal and opposite infall 
momenta, with a scale height of 200 pc. The global average 
dispersion from the mean plane will likely be rather larger, 
giving the first generation of stars an inflated final 
distribution. 

The critical test for the basic component of this model will 
be the observation at high redshift of 10** M d HIT regions 
emitting 10^ ergs per second in Lyman a, powered by 10^ ergs per 
second of 0 star radiation. They may otherwise resemble modern 
galaxies only in radial extent. Conversely, the most important 
area for theoretical extension is in halo evolution. 

Turning now to the larger range of possibilities in Figure 2, 
the predicted values of column density are proportional to Mpl/3 
t c ~^^ £ while radii are proportional to M p *^ t 
The ratio r/o increases systematically with r unless there is 
a strong correlation among the three perturbation parameters. The 



Since 


rotation velocity is proportional to M p *^t 

Mp is the parameter with by far the largest acceptable range, the 

rotation velocity essentially measures Mp. (Notice on Figure 2b 

that the constant v rot lines span only half an order of 

magnitude in M p .) For perturbation masses in excess of about 
12 8 

10 * , the model becomes inappropriate. If SO galaxies have 

characteristically high surface densities and rotation velocities, 
they might have arisen from these more massive perturbations. 

Similar models adequately extend this approach to low 
perturbation masses, providing either low surface density or very 
compact systems, depending on e. 

The largest problem the model has with low e perturbations is 
in deciding what happens next. Small systems, typically of high 
surface density, form and evolve rapidly in the core of the 

perturbation. What then happens to the rest of the perturbation, 

3 

with up to 10 times more mass? The result may depend on a fourth 
perturbation parameter having to do with the lurapiness of the 
original perturbation and the defocussing of the infall that 
lurapiness brings about. These less sheared systems could lead to 
elliptical galaxies with or without a small embedded disk, or the 
larger ones could become clusters of smaller ellipticals with a 
peculiar central object. Then again, they may turn out to lead to 
nothing ever seen, throwing some small doubt on either the model 
itself or on the regime of perturbation space available to the 



early universe* Once again it is halo evolution which is 
important, this time on a much grander scale* 

The reader will have noticed that the author has chosen to 
ignore the current enthusiasm for the idea that the apparently 
pervasive dark mass consists of matter in an unfamiliar form. The 
reason, if one is wanted, is that the hard' work to determine what 
sort of halo would evolve in the incredibly hostile environment of 
a self regulating galaxy has not yet been done. The dark mass may 
be in exotic forms, but then again it may not. More generally, 
one should be aware of the author’s innate prejudice against the 
recent notion that in evolving beyond the radiation dominated era 
the Universe became not only transparent but dull. 
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APPENDIX 


GROWTH, AMPLITUDE, SCALE, AND ANGULAR SIZE 
OF PROTOGALACTIC PERTURBATIONS 


A mass shell which reaches the origin at time t^ attains an 

2 

average minimum interior density P m ^ n = 37r/8Gt^ • A simple 
characterization of a perturbation is 


• <•*.«» [> - %) 2/3 ] 

so that the mass shell bounding Mp never quite falls in. The time 
dependence of infall is then 


t 0 2(t - t ) 

«M - Mp [1 - (^) 2 ] 3/2 - Mp ( J — ) 
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where the last form is that introduced previously for t - t « 

c 

For very early times in the development of the perturbation, 
the average density within a mass shell is 
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For the P^ n (M) dependence shown, substituting p •* 
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on 


the RHS above yields 
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Then substituting M ° y nr p a 4nr /(18nGt ) for M « Mp, the 
average density distribution close to the peak of the perturbation 
is 
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which has the desired quadratic dropoff with r that is 
characteristic of a rounded peak, justifying the P ra £ n (M) 
dependence chosen. The local density distribution is similar but 
with 1/3 replacing 1/5 in the quadratic terra. 

The relative amplitude of the perturbation is 

3 ( 3tt t -)2/3 _ 3 r P niin (0 \l/3 

5 1 2 t J 5 l p J 

c 

during the matter dominated phase. This can be written (for t « 

!!s , _J_ f„ f»i*/ 3 _i 2o_ 

p ^1/3 ^ t ' 1 + z 1 + z 

where fi = 2q Q , ty = H q and z is the redshift. The scale of the 

perturbations would be roughly the radius containing mass Mp at 
2 

density l/(6irGt ) 
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1 r 2M P Gt H ^1/3 3Mpc 

l n ) 




,1/3 


1 + Z 


1 + Z 


12 J 
10 M 


Since these perturbations are very close to our horizon, their 
presently observed angular diameter would be 


0 


r(l + z) n Q n 
c 


or roughly 10 (M p /10 12 M )^ 3 seconds of 


arc. 


independent of z. 
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FIGURE CAPTIONS 


Fig. I. The general properties of feedback controlled 
models of disk galaxy formation. (a) Infalling material 
converges toward the plane. (b) Owing to externally induced 
shear in the flow field, the convergence is to a disk. The disk 
material seeks its equation of state, phases, distribution, and 
star formation rate. (c) The disk grows in surface density and 
radius (higher angular raomentera falls in progressively later 
although continually accompanied by lower angular momentum 
material). The scale height shrinks with increased self 
gravity. Star formation and evolution contribute to an outflow of 
energy through the infalling material. At early times there is 
more material in the ’’near halo" (dashed outline) than in the 
disk, but the ratio decreases with time. (d) At this epoch, the 
mass in the near halo, AM has dropped to the point that it can be 
hydrostatically supported by the pressure associated with the 
outwardly diffusing energy. Infall continues at the outer edge 
but for the most part a quasi-static disk-halo system is 
established. It must, however, evolve dramatically as the disk 
exhausts its gas supply by star formation. (e) A halo which is 
shattered to stellar (and smaller) fragments at least as rapidly 
as the nuch denser disk will lead to the familiar system shown. 

a 

Such rapid fragmentation of the halo almost certainly requires 



inducement by energy leaving the disk* 

Fig. 2. The family of galaxies deriving from the HIT disk 
model. (a) The average disk column densities and radii are shown 
versus two composite parameters, t fi ow an< * ^t^. Both are 
constants for the evolution of a particular system. The former is 
a measure of the overall evolution rate, while the latter measures 
the rapidity with which the disk reaches Lyman a control. 

(b) Normalizing to a formation epoch at t^ = 10^ years, the family 
is shown versus the two remaining perturbation parameters. The 
perturbation mass is Mp, the shear parameter e. Primed values 
pertain to the normalization described in the text. Loci of 
constant surface density, radius, and rotation velocity are 
shown. The right hand axis indicates the fraction of the initial 
perturbation which enters the central galaxy by the time 
Lyman a stabilizes the halo. Inclusion of the near halo and edge 
growth raises the incorporation by roughly 3 to A. Notice the 
extreme differences in the breadths of ranges allowed the two 
parameters e and Mp* 
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